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We investigate differential electron momentum distributions in non-sequential double ionization 
(NSD1) with linearly polarized, few-cycle pulses, using a classical model based on a laser-assisted 
inelastic (e~ , 2e~ ) rescattering mechanism. We show that these yields, as functions of the momentum 
components parallel to the laser polarization, are highly asymmetric and strongly influenced by the 
absolute phase, i.e., the phase difference between the carrier oscillation of a few-cycle pulse and its 
envelope. Indeed, around a critical phase, such distributions change their sign in a radical fashion. 
This phase dependence provides a possibility for absolute-phase measurements which is, in principle, 
superior to the schemes involving high-order harmonic generation or above-threshold ionization. 



Recently, few-cycle laser pulses of intensities around 
or higher than 10 14 W/cm 2 have proven to be extremely 
important. A particular characteristic of such pulses is 
that they may have very high intensities and, still, carry 
much less energy than their longer counterparts, such 
that, effectively ionization is reduced. This has extended 
the damage threshold of solid-state materials up to the 
intensities in question , and has made the generation of 
high-order harmonic radiation up to astonishingly high 
frequencies possible 2]. Furthermore, their length, of the 
order of a few fs, permits controling processes such as 
molecular motion or chemical reactions , as well as the 
production of isolated, X- ray attosecond pulses Q| . 

In this pulse-length regime, the so-called "absolute 
phase", i.e., the phase difference between the envelope 
of the pulse and its carrier frequency, has considerable 
influence on several strong-field phenomena, as for in- 
stance high-order harmonic generation (HHG) [j| and 
above-threshold ionization [|| (ATI) 0, @ • Indeed, this 
phase determines features such as the maximal harmonic 
or photoelectron energy, the time profile of both phenom- 
ena, and the intensity of the ATI or HHG yields. 

This is not surprising, since the physics of HHG and 
ATI is directly related to the instantaneous, time depen- 
dent field. HHG, for instance, is described by a three-step 
mechanism in which an electron leaves an atom at an in- 
stant t{j through tunneling ionization, propagates in the 
continuum under the influence of the external laser field 
and, at a later time t\, recombines with a bound state of 
the parent ion, generating harmonics @. Slightly differ- 
ent processes, either in which elastic rescattering with the 
parent ion is taken as the third step, or in which the elec- 
tron reaches the detector without recolliding explain the 
high-order and low-order ATI peaks, respectively [Io| . 

In order to interprete the experimental data obtained 
in such cases, the precise knowledge of the absolute phase 
4> is required. This poses a serious practical problem, 
since this phase is difficult to stabilize, to control or to 
measure 111. For this reason, schemes for measuring <h 
have been suggested and realized, as, for instance, using 



the asymmetry in ATI photoelectron counts 0. 

In this Letter, we propose laser-assisted nonsequen- 
tial double ionization (NSDI) as a tool for absolute 
phase diagnosis. This phenomenon is being the sub- 
ject of very active discussion, which was triggered by 
differential measurements of electron momentum distri- 
butions performed with the COLTRIMS (Cold Target 
Recoil Ion Momentum Spectroscopy) technique, for lin- 
early polarized fields of intensities of the order of / ~ 
10 14 — 10 15 W/cm 2 incident in rare-gas samples [d|- Such 
measurements revealed very peculiar features, namely 
two symmetric peaks at p^u — p 2 \\ — ±2-jTIp, in the 
(Pi||>P2||) plane, where Pj\\(j — 1,2) and U p denote the 
momentum components parallel to the laser field polar- 
ization and the ponderomotive energy [l^ | , respectively. 

These features are explained by a physical mechanism 
very similar to those in HHG and high-order ATI. The 
main difference lies on the rescattering process at t\, 
which is now inelastic: the first electron gives part of 
its kinetic energy upon return to the second electron, so 
that it can overcome the ionization potential of the singly 
ionized atom and reach the continuum 

This laser- assisted rescattering process has been con- 
sidered b y se veral groups, using either classical [LJLsemi- 
classical jlol Uli Il7l ITq or quantum-mechanical |l9l I20L 
l2lj approaches, different types of electron-electron inter- 
action [jja, [jj], [lg , and neglecting or including electron- 
electron repulsion in the final states 0, 0, ^3 • So far, 
since the pulses involved were relatively long, they have 
been mostly approximated by monochromatic fields. 

Particularly what classical models concern, it is aston- 
ishing how well they reproduce the main features cither 
observed experimentally, or obtained by means of more 
refined, quantum mechanical methods. Indeed, recently, 
we have computed NSDI yields considering rescattering 
in its simplest form, i.e., electron-impact ionization, both 
classically and within a quantum-mechanical S-Matrix 
framework, with practically identical results [l7L ITsf . 

In this work, we use a similar classical model as in 
[rH Hs[ . in which an electron ensemble is subject to a 
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few-cycle pulse E(t) = —dA{t)/dt. The vector potential 
is given by 



The electron momentum distributions then read 



A(t) = A sin 2 (fM/2) s'm(ujt + ^)e a 



(1) 



with frequency to, amplitude Ao, absolute phase <f>, and 
£1 = Lu/n, where n denotes the number of cycles. The 
electrons are ejected in the continuum at a time to with 
vanishing drift velocities and from the origin of the coor- 
dinate system. The start times are uniformly distributed 
and the ejection probability per unit time, unless stated 
otherwise, is given by the quasi-static |22j tunneling rate 



R(t ) 



1 



\E(to) 



exp 



-2(2|£ 01 1) 3 / 2 
3\E(t )\ 



(2) 



where |£bi| is the ionization potential of the atom in 
question. Subsequently, these electrons propagate under 
the influence of only the laser held. Finally, some of them 
return to the site of their release and free a second en- 
semble of electrons through inelastic collisions at a later 
instant t\. 

The equations of motion of each pair in such electron 
ensembles, in atomic units, read 



[k + A(i )] 2 = 0, 



/ [k + A(t)} 2 dt = 

J t 



and 



Y^iPj + Mh)} 2 = [k + A(h)} 2 - 2\E 02 \ 



(3) 



(4) 



(5) 



Eq. © gives the energy conservation at the ionization 
time. Eq. J3J imposes restrictions upon the intermediate 
electron momentum k such that the electron returns to 
its parent ion. Finally, the third expression (Eq. (JSJ)) 
yields the energy conservation at the recollision time t\ . 
Thereby, the first electron gives part of its kinetic en- 
ergy E re t(ti) — [k + A(ti)] 2 /2 upon return to the second 
electron, so that it is able to overcome the ionization 
potential |_Eq2 | of the singly ionized atom. In terms of 
the momentum components parallel and perpendicular 
to the electric- field polarization, denoted p^-ii and Pj_l, 
respectively, Eq. (J5J is written as 

2 2 

5>i|| +A(i 1 )] 2 = [k + A(t 1 )] 2 -2| J B 2|-EpL- ( 6 ) 



Eq. 10 describes a circle in the (pi||,P2||) plane, centered 
at A(ti) and whose radius depends on E ret , \E 02 \ and on 
p 2 _i_(j = 1,2). The transverse momenta effectively shift 
the binding energy which must be overcome such that, 
depending on this quantity, there are situations for which 
the rescattering process is classically forbidden. 



dt R(t )S E xet (h) 
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where the argument of the S function gives the energy 
conservation upon return. The transverse momenta are 
integrated over. Details about this model are given in 



These distributions are displayed in the upper panels 
of Fig. 1, for neon j2^|, as density plots in the (pi||,P2||) 
plane. Their circular shapes and the maxima along 
Pill = P2\\ are features also present for monochromatic 
driving fields |24|, and mean, physically, that both elec- 
trons are leaving the parent ion most probably with equal 
parallel momenta. However, the fact that the yields are 
concentrated in only one quadrant of the (pi||,J>2||) plane, 
makes them strikingly different from the former distribu- 
tions, which are symmetric in (pi||,f>2||) ( — Pill , — P2\\)- 
Furthermore, for a narrow phase interval around a criti- 
cal phase 4> c (c.f. panel (c)), the sign of the momenta Pjn 
change in a rather radical fashion, and the whole yield 
is shifted from the first to the third quadrant. For in- 
creasing pulse length, these effects get less pronounced 
and practically disappear (Figs. l(e)-l(h)). The distri- 
butions then become symmetric and phase- independent. 

Important questions concern the physical origin of 
both the asymmetry and the critical phase: are they 
caused by the phase space or by the tunneling rate 
Depending on the parameters, a whole phase-space re- 
gion may become classically inacessible, such that the 
radius of the circle described by Eq. © would collapse 
and the corresponding NSDI yields would vanish. The 
quasi-static tunneling rate, on the other hand, favors the 
start times to for which the instantaneous field amplitude 
|-E(io)| is large, as compared to those for which |-E(£o)| 
is small. Thus, the contributions to the yield from the 
former or from the latter case would be enhanced or sup- 
pressed, respectively. 

In order to single out the influence of the phase space, 
we assume that the electrons belonging to the first en- 
semble reach the continuum at a constant rate. Such 
results are shown in panels (i)-(l) of Fig. 1, and are 
radically different from those obtained with the more 
realistic, quasi-static tunneling rate. Indeed, the mo- 
mentum distributions, though asymmetric, exhibit peaks 
in both first and third quadrants of the (pi||,P2||) plane, 
vaguely resembling those obtained with monochromatic 
driving fields. The asymmetry is expected, since, for 
such pulses, the relation A(t) = —A(t ± T/2), and thus 
|r(ti,to,Pi||,J>2||)| - \T(t 1 ±T/2,t ±T/2,-p lh -p 2ll )\, 
where T = 2tt/uj, which was true for monochromatic 
fields, no longer holds. However, the huge effects ob- 
served in the upper panels are absent. Physically, this 
means that there is a momentum region for which the 
rescattering process is classically allowed but for which 
the probability that the first electron reaches the contin- 
uum is very small. Consequently, even if this region is 



FIG. 1: Electron momentum distributions along the laser po- 
larization, for neon (|-Em| = 0.79 a.u. and \Eq2\ = 1.51 a.u.), 
subject to pulses of intensity I = 4.7x 10 14 W/cm 2 and carrier 
frequency u — 0.057 a.u., respectively. In panels (a) to (d) 
and (i) to (1) we consider a four-cycle pulse (n = 4), whereas 
in panels (e) to (h) the pulse length is varied. Panels (a) to 
(h) and (i) to (1) were computed with the quasi-static and 
a constant tunneling rate, respectively. Panels (a) and (i): 
4> — O.ln; panels (b) and (j): <f> = 0.57r; panels (c) and (k): 
4> = 0.87r; panels (d) and (1): (f> — 0.9ir. In panels (e), (f), and 
(g) n = 8, n = 12 and n = 16 respectively, and (f> — O.ln. In 
panel (h), n is the same as in (g) and (f> = 0.9-7T. 
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FIG. 2: Parallel electron momenta py along pi« = p 2 ||, for 
transverse momenta pix = P2X = 0, as functions of the start 
times to of the electrons belonging to the first ensemble, to- 
gether with the quasi-static tunneling rates. The remaining 
parameters are the same as in panels (a)-(d) of Fig. 1. The 
times to aregiven in units of the field cycle. 



large, its contributions to the yield will be negligible. 

In Fig. 2 we analyze this effect in detail. Therein, 
the electron start times to are plotted, together with the 
quasi-static rate. We restrict the parameter range so that 
the classically allowed region is most extense |25| , taking 
parallel momenta along the diagonal || = P2\\ = P\\ and 
vanishing transverse momenta |16|. We consider only 
pairs (ti , to) of start and return times for the first electron 
such that its excursion time At = t\—to in the continuum 
is of the order of T/2 [2(|. For the specific parameters 
of the figure, there exist mainly two sets of electrons for 



FIG. 3: Parallel electron momentum distributions for pulses 
of intensities I = 4 x 10 14 W/cm 2 and / = 8 x 10 14 W/cm 2 
(upper and lower panels, respectively). The electrons were 
ejected with the quasi-static tunneling rate. The remaining 
parameters are the same as in panels (a)-(d) of Fig. 1. 



which the quasi-static rate is large and, therefore, whose 
contributions are relevant: those ejected at 1.5T < to < 
2T, with positive momenta, and those released at IT < 
to < 2.5T, with negative momenta |27j . 

For a large range of absolute phases <fi < (f> c , electron- 
impact ionization from the latter set of trajectories is 
classically forbidden. Thus, the distributions concentrate 
on the first quadrant. Around the critical phase, this pro- 
cess is allowed for both sets of electrons and the tunneling 
rates are comparable. Consequently, there are relevant 
contributions to the yield in the first and third quadrants. 
This changes for larger phases, as for instance cf> = 0.9-7T. 
In this case, the electrons ejected at 2T < to < 2.5T are 
favored and Pj\\ are mainly negative. 

The role of the phase space is even more evident as the 
driving-field intensity is varied, as shown in Fig. 3. By 
doing so, one is changing the radius of the circle described 
by Eq. ©, and thus the region in the (pi||,P2||) plane for 
which rescattering is classically allowed. Therefore, the 
critical phase may change. 

For a lower intensity than that in Fig. 1, the yield 
in the negative momentum region appears for a phase 
larger than (j> c — 0.8ir (c.f. panels (c) and (d)). This 
is due to the fact that the classically allowed region for 
pu < is almost vanishing. Thus, even if the tunnel- 
ing rates are comparable, the first electron, upon return, 
no longer possesses enough kinetic energy to release the 
second electron in a way that both leave with negative 
parallel momenta. For a higher intensity, apart from the 
fact that both regions are classically allowed, the proba- 
bility that the second electron is released with negative 
parallel momentum is larger. Therefore, the transition 
occurs where expected, as displayed in panel (g). This 
is confirmed by Fig. 4, where, as the intensity decreases, 
the dominant set of ionization times (2T < to < 2.5T), 
corresponding to p|| < 0, collapses. 

In conclusion, we perform a theoretical investigation 
of NSDI with few-cycle pulses, using a classical model 
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FIG. 4: Parallel electron momenta py along — p 2 || and for 
vanishing transverse momenta, as functions of the start times 
to, together with the quasi-static tunneling rates, for several 
driving- field intensities and the critical phase <j> — 0.8ir. The 
remaining parameters are the same as in Fig. 3. For the 
lowest intensity, rescattering caused by electrons ejected at 
2T < to < 2.5T is classically forbidden, so that the corre- 
sponding curve is not displayed. 



based on electron impact ionization. Both electrons have 
equal final momentum components parallel to the field 
polarization, which are mainly positive or negative. The 
sign of such momenta and their most probable values 
depend on the absolute phase <j). In particular, around a 
critical phase, these momenta change sign. Such features 
are explained as the interplay between the tunneling rate 
for the first electron and the phase space. 

The changes in the NSDI yields upon a critical phase 
are far more extreme effects than those observed for 
HHG or ATI. In fact, nonsequential double ionization 
has an advantage over the other two phenomena: the 
phase space region contributing to the process is con- 
fined. Thus, for particular ranges of <p, it is easier to 
make a whole region either classically forbidden, by mak- 
ing the radius in Eq. (6) collapse, or irrelevant, by reduc- 
ing the corresponding ionization rate. Therefore, NSDI 
is, in principle, a tool for absolute phase diagnosis which 
is more efficient than the existing schemes. 

We thank A. Fring for useful discussions. This work 
was financed in part by the Deutsche Forschungsgemein- 
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